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Abstract 


The influence of systematic errors in the parameters of the models describing 
the geometry and the atmosphere on the profiles of trace gases retrieved from simulated 
solar occupation spectra, collected at satellite altitudes, is investigated. It is concluded 
that, because of smearing effects and other uncertainties, it may be preferable to 
calibrate the spectra internally by measuring absorption lines of an atmospheric gas such 

as CO^ whose vertical distribution can be assumed rather than to rely on externally 
supplied information. 
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Introduction 

It i3 expected that high-resolution, infrared, solar, occupation spectra covering 

wide spectral regions will soon he collected from satellite altitudes/ These spectra will 

add significantly to our knowledge of atmospheric properties, including the latitudinal and 

vertical distributions of many trace gases. However, the quantitative analysis of these 

spectra requires a knowledge of the geometry of- the experiment and of the 

characteristics of the atmosphere along the ray paths. Although this information can be 

obtained from external sources, features in the collected spectra, including, for example, 

those of gases such as CO , whose vertical mixing-ratio profiles are well known, can be 

2 

2 

used to calibrate the spectra Internally. 

In this paper the influence of systematic errors in the path geometry and the 

atmospheric characteristics on the derived mixing-ratio profiles of trace gases is 

2 

estimated with the simulation methods used previously. It is assumed that the data 

consist of the equivalent widths of weak lines. The relative merits of internal and 

external calibration of spectra obtained from satellite altitudes are compared and it is 

concluded that, in many cases, it may be preferable to calibrate the spectra internally by 

measuring lines of a gas such as CO^ whose vertical and latitudinal distribution is 

4 2 

reasonably constant. Fanner et al. have previously found that calibrations, based on 

CO lines, of solar spectra obtained from balloon altitudes also reduced systematic errors 
2 

associated with the atmospheric models and observational geometry. 
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Systematic Errors Associated with the External Calibration Method. 

The equivalent width of a weak atmospheric line can he written as 


W( Zl ) = E SjPjdjSl^ 


(1) 


where S is the line intensity, p is the mixing-ratio of the absorbing gas, d is the air 

} } i 


density, and 51 is the path length in the jth atmospheric layer traversed by the ith ray 
ij 


which reaches the instrument at an angle z^ with respect to the local zenith. Some of the 


quantities affecting these equivalent widths are summarized in Table 1. In general, if the 


line intensity depends on parameters s , the mixing ratio on parameters t , the density on 

g p 

parameters x , and the incremental path on parameters y , then the change in equivalent 
q. r 

width caused by changes in these parameters 


dW(^) = £ Sj pjdj 61 ^ 


as. 
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It is convenient to assume, initially, that the line intensity is known accurately 


and is independent of temperature. In this case the amount of absorber a(z^) along the 


incoming rays 


a(z ) = Wfe )/S , 
i i 


(3) 


and 
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where we have assumed that the set of x variables is independent of the set of y 

q r 

variables and the set of t variables is the union of the x and y sets. 

<1 q t 


The effects of errors in d and 61 on the values of p retrieved from a set of 

i ij i 


afz 5 values have been studied by creating a set of measurements a(z ) with an appropriate 
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set of values for the parameters x , and y , Mixing ratio profiles were retrieved 

q r 

successively from these data using incorrect values for each of the parameters x or y . 

q r 

If, for example, the incorrect value x^ + Ax , was used then, since the correct set of 
absorber amount values were analyzed, A a (z) - 0, and Bq. (4) becomes 


E 

This can be written as 


life. 


■I 


3x, 


54 


Ax, 


= 0. 


e f&* fi =0. 

i 0, d , 

If the local mixing-ratio values are independent then 


fa *fi ... 

p, 


similarly 


+ = 0. 
<5. 


(5) 


J 


i) 


5,6 


The program described by Snider and Goldman was used to calculate the 
absorber amounts along the paths of 125 rays, with nearly equally-spaced tangent heights 
between 0 and 50 km, to a satellite 500 km above the surface and at 45° latitude. The 


absorber was assumed to have a constant mixing-ratio in all three dimensions and to be in 

an atmosphere with the temperature and density profile of the G.S. Standard 
7 

Atmosphere. A wavelength of 5 pm was assumed. 

For these conditions the distance along; a. ray, such as that shown in Fig. 1, from 
the tangent point to the satellite is about 2500 km and the distance to a point 50 km 
above the tangent point is about 800 km. Variations in the effective radius of the earth or 
in the atmospheric temperature, density, and mixing-ratio profiles along these paths have 


3 
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not been modelled, nor have the effects of the finite FOV of the instrument and satellite 

6 

motion during the period of data collection. Kantor has shown the variability of the 
density of the atmosphere between 20 and 60 km above the surface from the equator to 
65° over periods of up to 72 hours and distances of up to 370 km is typically several 
percent. The finite angular FOV causes a spread of the beam at the tangent point of 
about 2.S kra/mrad. Satellite motion may cause the tangent heights of the rays to change 
at rates as high as 3 km/sec. Thus the vertical resolution of typical measurements made 
with a FOV of 1 mrad and for observation times of about one second is limited to a few 
kilometers. These effects also set upper limits to the accuracy of the information which 
can be derived from the spectra. 

Layer Thickness 

The paths of rays through the atmosphere and the total air mass traversed were 
calculated by dividing the atmosphere into layers with the same vertical thickness. The 
dependence of these calculated air masses on layer thickness is shown in Fig. 2 for 
thicknesses of from 25 m to 20 km and rays with tangent heights of 2.6, 12, and 48 km. 
The calculations are made by assuming there is a constant vertical temperature gradient 
in each layer. For the two lowest rays, this is a good approximation for all layer 
thicknesses. The larger variations in the calculated air mass traversed by the ray with a 
tangent height of 48 km are attributed to the rapid changes In the temperature gradient 
of the Standard Atmosphere near 50 km. In the rest of this work a layer thickness of 0.25 
km v/as used. This still caused irregularities seen in some of the figures near the 
tropopause level. 


4 


original page 18 
OF POOR QUALITY 


Atmospheric Refraction 

The atmospheric dispersion from the visible to beyond 20 pm is small and, to a 

very good approximation, rays of all wavelengths can be assumed to travel the same 

paths. However, if refraction is not considered, the calculated air masses traversed by 

rays with tangent heights between 0 and 90 km have the errors shown in Fig. 3. If the 

mixing ratio profile is retrieved by neglecting refraction, the calculated path segments 

61 required in Eqs. 1 and 5 are incorrect and errors in the retrieved mixing-ratio profiles 
Ij 

are produced, as shown in Fig. 4. The irregularity near 10 km is due to the rapid change in 
the atmospheric temperature gradient at the tropopause. 


Reference Pressure Errors 

The atmospheric pressure and density profiles were obtained from the vertical 

temperature profile and the sea level pressure P , If the correct temperature profile is 

o 

used, the pressure and density at any other level and the air mass traversed by a ray with 

a fixed tangent height are proportional to P . Thus 

o 

Ad/d = AP /P , (6) 

o o 

and hence from Bq. (5) 

Ap/p = -AP /P . (7) 

o o 

The fractional errors in the retrieved mixing-ratio, A p/p, shown in Fig. 5, and 
obtained for the case where the ray path geometry is known accurately but AP /P is 

O G 

equal to ±5%, are not independent 6f height, as expected from Eq. (7), because the change 
in the air density causes the refraction effects to be modelled incorrectly. 
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Incorrect reference level pressures may be used because of uncertainties in the 

large-scale, global variations in atmospheric pressure, density, and temperature. Xn 

addition, smaller-scale variations in these quantities of several percent may occur over 

8 

regions much smaller than those traversed by the rays. These variations are difficult to 
model in the spectral analyses described here and they may limit the accuracy of the 
retrieved mixing-ratio profiles. Since pressure and height are related, the use of an 
incorrect reference pressure at a known altitude is equivalent to choosing an incorrect 
altitude for the reference pressure. These relations can be estimated for an isothermal 
atmosphere in which the pressure and density are given by 


P(h) = P exp (-h/H), 08) 

o 

and dOh) = d exp C-h/H), (9) 

o 

where the scale height of the atmosphere 


H = RT/Mg, 


( 10 ) 


is approximately 7 km below 50 km in the Standard Atmosphere. 
Since 


Ah/H = AP /P , Oil) 

o o 

a five percent error in P is equivalent to an error of 0.4 km in the tangent heights of all 

o 

the rays. 


Errors in Atmospheric Temperature 

The dependence of atmospheric density on temperature can be approximated by 

considering an isothermal atmosphere with a fixed reference pressure P . Prom ,Eqs. (8) 

o 

and CIO) 


6 


AP = MghPAT/RT 2 , 
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( 12 ) 


and* 

for a perfect gas, 



d » MP/RT. 


Thus 

2 

Ad = MAP/RT - MPAT/RT , 

(13) 

and 

Ad/d = (h - H)AT/HT. 

(14) 


As the temperature of the entire atmosphere is changed the density remains 
constant near h = H hut it increases at a rate of about 2.596/K near 50 km. Prom Eq. (4) 
the dependence of the errrors in the retrieved mixing-ratio profiles on these density 
changes is also described by Sq. (14). 

The fractional changes in the atmospheric density caused by changing the 
temperatures at all levels in the standard atmosphere by ±5%K and the corresponding 
fractional errors in the retrieved mixing-ratios in Fig. 6 are similar to those predicted by 
Bq. (14). Larger errors in the mixing-ratio may occur if the line intensities depend on 
temperature. It is planned to discuss these effects elsewhere since they may be used to 
obtain improved estimates of the temperature profile if the mixing-ratio of the absorbing 
gas is known. 

The differences in the tangent heights of rays which pass through the same air 
masses in isothermal atmospheres at different temperatures are obtained from Bqs. (9) 
and (14) 

Ah = (H - h)AT/T . (15) 

These differenc.es increase to about 0.4 km near 50 km for a one percent change 
in the temperature of the entire atmosphere. In addition to mixing-ratio errors caused by 

i " 
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the use of incorrect vertical temperature profiles unmodelled horizontal temperature 
variations may also affect the accuracy of the retrieved mixing-ratio profiles. 


Satellite Heights 

It is seen from Fig. 1 that a change Aa in the satellite height causes a change 

Aa cos© — Aa in the tangent heights of all the rays. From Bq.. (9) Ad/d + Ah/H = 0. 

Thus, for an isothermal atmosphere, the rate of change of density with height, and hence 

of the air mass traversed with the tangent height of a ray, is about 14%/km. The errors in 

the retrieved mixing- ratios caused by errors of ±1 km in the satellite height, shown in 

Fig. 7, are approximately 14% and the vertical dependence is very similar to that caused 

by errors in P . To reduce the errors to less than 1% the satellite height must he known 
o 

to better than 100 m. 

i 

Latitude Error 

r 

The effective radius of the earth and the acceleration due to gravity are both 
latitude dependent and the air masses traversed by rays reaching an observer at fixed 
incoming angles vary with latitude. A series of calculations has shown that, for an 
observer at SO km, the air masses traversed by the rays varied by less than ±4% from the 
value at 45°, as the latitude was varied from 0° to 90°, and the changes were essentially 
the same for all rays with tangent heights between 0 and 50 km. Thus, since the latitude 
is always well determined,' these changes should not be a significant source of error. . • 
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Incoming Angle Errors 

Atmosphere paths are calculated from the angle z shown in Fig. 1, between the 

local zenith and the arriving rays. This angle can be measured directly or it can be 

calculated from the known position of the instrument and the time of observation, 

provided atmospheric refraction is included. 

An error Az in this angle causes an error Ab * LAz in the tangent height, where 

L is the satellite-to-tangent height distance. For a satellite at 5D0 km, Ab/Az is 

approximately 2.4 km/mr or 40 km/l°. This corresponds to a rate of change in the 

atmopheric density at the tangent height of about 600%/1'\ The changes in the retrieved 

mixing-ratio caused by systematic errors in z of 0.01° and corresponding errors in b of 
* 

0.4 km are shown in Fig. 8. As expected, they are nearly independent of height and 
average about 6%. 

When the angle z is calculated from the parameters of the satellite orbit, errors 
may arise from incorrect estimates of the satellite height, the position in its orbit, the 
plane of the orbit, or the time of observation. For the two-dimensional geometry in Fig. 1 
z = 90 + 0, (16) 

Az = A0 . 

cos© o (R + b)/(R 4 a) , (17) 

Ab = LAz , 

2 2 

and -sin0A0 - -(R + b)Aa/(R + a) + (a - b)AR/(R 4 a) + Ab/(R 4 a). 

For a satellite with a = 500 km, R = 6400 km, and b « a, errors of 0.5 km in 
a, of 6 km in R, or of 0.01° in 0 or z cause 7% errors at all heights in the derived 
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mixing-ratio profile. The same error is caused if the position of the satellite in its orbit 
is uncertain by 1.2 km or if there is a timing error of 0.16 sec. Since all of these errors 
affect the derived mixing-ratio profiles in the same manner as errors in the reference 
pressure* additional information must be supplied if the origins of suspected errors are to 
be identified. 

Satellite Motion 

For the two-dimensional geometry in Fig. 1 and for a satellite in a circular orbit 

at 500 km the tangent height of a ray changes at a rate of about 2.8 km/sec, in the 

absence of refraction. The distortion of these paths due to changes in the atmospheric 

9 10 

density gradient and to refraction have been discussed by Garriott and Russell. 

The corresponding rate of change of air mass* in an isothermal atmosphere* 
obtained from Eq. 9 is about 40%/sec. These rates of change of air mass as functions of 
the change in apparent zenith angle were also calculated for an atmosphere with a 
temperature gradient of 1 K/km and for the U. S. Standard Atmosphere with refraction 
effects included. The results are shown in Fig. 9. In both cases the mean rate is about 
40%/sec. However the rate varies significantly with the tangent heights of the rays 
because it is affected both by refraction and by the atmospheric temperature gradient 
near the tangent height. Since the mean air mass traversed by the sets of rays reaching 
an instrument is obtained by weighting over the POV these results show that the weighting 
is a function of the temperature profile and changes in the temperature gradient must be 
considered when the mixing-ratio profile is obtained from actual data. 
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The amount of an absorbing gas along the atmospheric paths traversed by solar 
rays can be obtained from the equivalent widths of weak absorption lines in solar spectra 
provided the models used to describe the geometry of the experiment and the structure of 
the atmosphere are known. The effects of systematic errors in these parameters on the 
retrieved mixing- ratio profiles of the absorber have been estimated. With the exception 

V 

of the effects of atmospheric refraction and errors in the temperature profile these errors 
produce nearly identical effects on the retrieved profiles. These effects are equivalent to 
an error in the tangent heights of the rays producing the set of spectra which is 
independent of the height of the ray. The magnitudes of these effects agree with those 
estimated by making simplifying assumptions about the experiment design. 

The existence of these errors can be detected by comparing the mixing-ratio of 
a gas whose vertical distribution in the atmosphere is known with that obtained by 
analyses of the equivalent widths of its lines in the spectra. However, because of the 
similarities of these effects it is not possible to determine, from these profile differences, 
either the number of parameters in Table 2 which are in error or the magnitudes of these 
errors unless additional information is supplied or different types of measurements are 
made. No discussion of the influence of several smearing effects present in the spectra 
on the derived mfcdrig-ratios has been given although it has been pointed out that these 
can substantially reduce the vertical resolution. 11 Incorrectly modelled smearing 

effects may lead to systematic errors substantially larger than those in Table 2. In the 

* 

remainder of this discussion the method of Internal calibration is discussed. Although 
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itself subject to systematic errors these are expected to be smaii. This method does avoid 
many of the errors due to the incorrect modelling of smearing effects as well as many of 
the systematic errors discussed above. 
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Internal Calibration of Occultation Spectra 

When the information about the atmospheric paths is obtained from external 

sources several errors affect the retrieved profiles in nearly identical ways. This occurs, 

in part, be ause the spectral information concerning the atmosphere, carried by a ray 

with a fixed tangent height and which originates outside the atmosphere, and is collected 

outside the atmosphere, Is not influenced by the positions of either the source or the 

detector along the ray and hence by the particular sets of parameters which describe 

these positions. It has also been noted that the quality of the spectra is degraded by the 

11 

finite FOV of the instrument, satellite motion, atmospheric refraction, and temperature 

12 

- and density - profile variations along the atmospheric path and that it is unlikely that 
all of these effects can be satisfactorily modelled in the analysis. 

Internal calibration removes many of these uncertainties but also introduces new 


ones. In this method, information about the atmospheric path is obtained from absorption 

lines in the spectra. Lines of CO are often measured because its vertical and latitudinal 

2 

mixing-ratio is nearly constant and the parameters of many CO lines are accurately 

2 

known. We consider, initially, the information obtained by measuring the equivalent width 


of a single weak CO line, with a temperature-independent line intensity, in a single 
2 

spectrum. This allows one aspect of the path, such as the total amount of CO , to 

2 
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be obtained. The additional information obtained by measuring the equivalent widths of 
several lines in different regions of the curve of growth, by measuring lines whose 

t 

intensities are temperature dependent, or by measuring the shapes of lines is considered 
later. 


If W and W are the equivalent widths, in a single spectrum, of weak lines of 
0 1 

CO and a trace gas, respectively, and if the corresponding line intensities S and S are 
2 0 1 

temperature independent, then, from Eq. (3), the amount of CO along the path 

2 


a = W / S , 

0 0 0 

and the amount of the other absorber 


(18) 


a = W /S 
1 11 


(19) 


If the CO mixing-ratio p is constant with height and along the path traversed 
2 0 

then the total air mass traversed 


A = W/Sp , (20) 

0 0 0 

and the mean mixing-ratio of the absorber along the ray path, relative to that of CO 

2 

p = a /A . (21) 

11 

These simple results are obtained because the incremental equivalent widths are additive. 

hi this approach both equivalent widths have been averaged in some (unknown) manner 

over the FOV of the instrument and they are weighted by the other factors discussed 

above which affect the quality of the data. If the spectra are obtained by a Fourier 

Transform Spectrometer information concerning both of the lines W and W is collected 

0 1 

simultaneously. Provided the lines are both weak and neither of the mixing-ratios change 

r * i •' 1 ' 

rapidly over the FOIT both of their equivalent widths axe subjected to the same weighting 
processes which are then automatically ratioed out cjf the final results. Thus the values 
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for the total air mass traversed and the mean mixing-ratio of the unknown absorber given 

by Eqs. (19) and (21) do not contain the implicit uncertainties associated with the 

uncorrectable smearing effects of the external calibration methods. 

Additional measurements or assumptions are still required to obtain the 

3 

distribution of the unknown absorber with height, in the previous paper this distribution 
was obtained from the measured equivalent widths and the known geometry and the 
atmospheric structure. Thus the mixing-ratios could be determined either as a function 
of altitude or of pressure. 

in the internal calibration method the geometry is inferred by comparing the W 

0 

measurements with another set obtained from Eq. (20) and by assuming the atmospheric 

temperature, density, and CO profiles are known. The unknown mixing-ratio is then 

2 

3 

obtained as described previously. These atmospheric paths will differ from the actual 
paths if these assumed profiles are incorrect. Thus errors in the mixing-ratio profile of 
the unknown absorber will occur which will depend on the nature of the systematic error 
and the method used to describe the profile. For example, an error of APo (e.g. of one 
percent) in the reference pressure Po does not affect the estimates of the total air mass 
traversed since air masses are derived from Bq. (20). However, from Eq. (11) it can cause 
a corresponding error Ah ( — 0.1 km) in the tangent heights of all rays and hence the entire 
mixing-ratio is displaced vertically by Ah. 

The effects of other types of errors can be estimated in a similar fashion and 
they are listed in Table 3. With the exception of errors in the temperature profile these 
effects are independent of the tangent heights of the rays. If the values in Table 3 are 

r ' I ■ * 
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added in quadrature the total uncertainty is less than 1 km. This is less than the vertical 

resolution due to the smearing effects already discussed. It is also smaller than the 

vertical resolution of the trace gas profiles expected to be obtained from mid-infrared, 

occultation spectra obtained with state of the art Fourier transform spectrometers at 

balloon or satellite altitudes. If reasonable assumptions are made concerning the rate of 

data collection and the SMS this resolution is typically several kilometers. Thus internal 

calibration is appropriate for these measurements. 

The internal calibration method does require the CO line intensities and its 

2 

mixing-ratio. Although the spectrum of CO has been extensively investigated the wide 

2 

range of intensities required to obtain accurate equivalent width information may require 

additional quantitative laboratory measurements of CO spectra. A careful evaluation of 

2 

4 

the global and temporal variations in the atmospheric CO mixing-ratio is also needed. 

2 

13 

Mauersberger and Finstead have summarized some of the information on atmospheric 

CO and have described some higher- than-expected CO mixing-ratios obtained above 
2 2 

25 km at night. 

Some additional types of information which may improve the accuracy of the 
results is discussed below. 

Relative Times of Data Collection 

The rate of change dz^/dt of the zenith angles of the incoming rays and the 
corresponding rate of change db/dt of the tangent height of rays reaching a satellite are 

nearly constant in a given occultation. Thus the scatter in Hie equivalent tangent heights 
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b^ with time allows the consistency of the data to be estimated. Nearly constant 
systematic shifts between all of the values obtained by internal calibration of the spectra 

and those calculated from the orbital parameters may be due, for example, to timing 
errors. Differences at low tangent heights below 10 km may be caused by improperly 

corrected refraction effect and uncertainties in the CO^ mixing-ratio or departures from 
LTE may produce differences between the observed and calculated curves for tangent 

heights above 70-80 km. However, good agreement can usually be expected for rays with 
tangent heights between about 20 and 50 km and these data can be used to improve the 
calibration of spectra which do not contain CO lines suitable for calibration. 


Measurements of Several Lines 

The bands and lines of CO^ have a wide range of intensities. Thus by measuring 
the equivalent widths of a number of lines in each spectrum the consistency of the 

measurements can be estimated and those lines which are too weak to be measured 
accurately or which lie outside the linear region, can be rejected. The effective SNR is 
also improved by measuring a number of lines. 

Because of the very large variation in the air masses traversed by rays with 

tangent heights from 0 to above 50 km it is necessary to measure CO^ lines with a 
corresponding range of temperature-independent intensities in order to maintain high 

precision. In general, this will require measurements of CO lines, with small lower state 

2 

energies, in bands with different intensities. 
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Lines Not in the Linear Region 

It has so far been assumed that the equivalent widths of the CO calibration 

2 

lines can be obtained from Eg. Cl), which can be written as 

W (z. ) = Z AW„ . 

1 £ ij 

If the CO lines are not in the linear region of the curve of growth the 
2 

incremental widths AW^ depend on the local atmospheric temperature and density and 
they are no longer additive. The weighting of such lines over the FOV and other smearing 

effects is different from that for weaker CO lines or weak lines of trace gases. 

3 

In these cases one of the principal advantages of the internal calibration method, 

the ability to assume that the smearing effects are the same for all lines, is lost. It is 

therefore important to choose CO calibration lines which are in or close to the linear 

2 

region of the curve of growth. 


Measurements of Temperature Dependent Lines 

Because the atmosphere is not isothermal the weighting of the equivalent widths 

of the lines over the FOV depends on the temperature dependence of the lines in addition 

to other effects previously considered. Since the strongest lines of mo3t trace gases are 

not strongly temperature dependent CO lines with a similar temperature dependence 

2 

should be chosen for calibration. Discussions of the atmospheric information contained in 

14,15 

temperature dependent CO lines are given elsewhere. 

2 

Curve Fitting 

Spectral curve-fitting is a powerful method of extracting information 
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If the models describing the spectrum are accurate. Park has shown that the 

effects of satellite motion during data collection and also other effects distort the 

occultatlon interferograms and consequently the spectra derived from them and that 

these smearing effects on the line shapes are difficult to model accurately. These 

3 

methods are not explored in this paper. Shaffer and Shaw have shown that more 
information on the vertical distribution of gases can he obtained by collecting as may j 
spectra as possible during an occultatlon than by inferring the distribution from the 
profiles of the lines in a single spectrum. One of the principal effects of measuring 
equivalent widths is to smooth the noise in the data analysed. 

Conclusions 

Internal calibration of satellite occultatlon spectra from measurements of the 

equivalent widths of weak, temperature-independent CO lines may he preferable to 

2 

calibration based on error-prone external information. It is expected that internal 
calibration should allow the mean tangent heights of the rays to estimated to about 1 km 

a 

and the mixing ratios of trace gases to about five percent. 

These uncertainties increase if the intensities of the lines used for calibration 
are temperature dependent and if tfye lines are not in the linear region pf the curve of 

I 

growth because the weighting of the absorber over the BOV and other smearing effects 
depend on the line parameters. 

However, because of the large .numbers and the wide range of intensities of CG 

2 

lines appearing in the solar spectrum the accuracy of the internal calibration method can 
be maintained over wide height intervals. 
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Table 1. Quantities Affecting the Equivalent 
Widths of Weak Lines in Solar Spectra 




HR 


Table 2, Sources of Error and Their Effects 

* 

on Retrieved Mixing Ratio Profiles 
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Quantity in Error 

Mixing Ratio Error, 
Dependence on Height 

Change Required for 
1% Error in Mixing Ratio 

Refraction 

only important below 25 km 


Reference Pressure 

independent of height 

AP » 10 mb 
o 

Temperature Profile 

increases with height 


Satellite Height 

independent of height 

Ah = 80 m 

Latitude 

ir ti it 


Radius of Earth 

II VI It 

1 km (at 45°) 

Angle of Incoming Rays 

it « ii 

0.002° 

• 

Time of obs. 

it n n 

0.03 secs 

Position in orbit 

I ii it it 

0.25 km 
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Table 3. Dependence of the Errors in the Tangent Heights of Rays 

* 

Made by Using Internal Calibration Methods on 
Systematic Errors in the Assumed Parameters Values 


Parameter 


Systematic 
Fractional 
Errors (%) 


Corresponding Change in 
Tangent Height of Rays 


Sea level pressure 


0.1 km 


Temp Profile 


1 (3K) 


<0.5 km 


0.1 km 


0.1 km 







C0 2 mixing ratio 
Line intensity 
Equivalent width 


1 


0.1 km 














Fig. 1 

Fig. 2 
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t 

Fig. 3 

Fig. 4 
Fig. 5 
Fig. 6 

Pi|7 
Fig. 8 
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Occultation Geometry. For a mean earth radius R = 6360 1cm and a satellite 

9 

height a = 500 1cm the distance L is about 2500 1cm, the angle 6 is about 21° and 

4 ( * 

the zenith angle z is about 111°. 

The dependence of the air mass traversed tgr rays with tangent heights of 
approximately 3, 12, and 48 km on the thickness of the atmospheric layers used 

i 

in the calculation. The fractional change in air mass with respect to that 
calculated for 25 m thick layers is shown. 

The dependence of the fractional air mass error along the rays, due to the 

Ir 

neglect of atmospheric refraction, on the tangent height of the rays for the 
case of an observer at 500 km. 

The height dependence of the fractional errors in the mixing-ratio profile 
caused by the neglect of atmospheric refraction in the retrieval. . 

The height dependence of the fractional errors in the mixing ratio profile 
caused by five percent errors in sea level pressure. 

The height dependence of the fractional errors in the retrieved mixing-ratio 
profile caused by assuming ±5 £ changes at all heights in the IT. S. Standard 
Atmosphere temperature profile. 

The height dependence of the fractional errors in the retrieved mixing-ratio > 
profile caused by errors of ±1 km in the satellite height. { >; 

i 

The height dependence of the fractional errors in the retrieved mixing-ratio 
profile caused by errors of ±0.01° in the angles of the incoming rays. b 
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Fig. 9 Fractional rate of change of air mass (sec ) for an observer at 000 km 

a) TX.S. Standard Atmosphere hi) atmosphere with a temperature gradient of 
1 K/km. The tangent heights of the rays decrease from 90 to 0 km. 
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